Glass transition line in Ceo: a mode-coupling/molecular-dynamics study 
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We report a study of the mode-coupling theory (MCT) glass transition line for the Girifalco model 
of C60 fullerene. The equilibrium static structure factor of the model, the only required input for the 
MCT calculations, is provided by molecular dynamics simulations. The glass transition line develops 
inside the metastable liquid-solid coexistence region and extends down in temperature, terminating 
on the liquid side of the metastable portion of the liquid-vapor binodal. The vitrification locus 
does not show re-entrant behavior. A comparison with previous computer simulation estimates 
of the location of the glass line suggests that the theory accurately reproduces the shape of the 
arrest line in the density-temperature plane. The theoretical HNC and MHNC structure factors 
(and consequently the corresponding MCT glass line) compare well with the numerical counterpart. 
These evidences confirm the conclusion drawn in previous works about the existence of a glassy 
phase for the fullerene model at issue. 



I. INTRODUCTION 

The onset of a glassy phase characterized by positional 
disorder in the Girifalco central potential model of Ceo [l| 
has been recently documented by some of us via Molec- 
ular Dynamics (MD) studies 0, Q. Interest in the vitri- 
fication process in fullerenes stems not only from the in- 
trinsic relevance of this class of materials, but also from 
the nature of their interparticle interaction. The Giri- 
falco potential appears in fact to be "marginally" short- 
range, giving rise to peculiar effects when we consider 
the interplay between the decay of the interactions and 
both the characteristics of the phase portrait and the 
glass forming ability of this model. In particular, follow- 
ing an initial debate on the existence of a stable liquid 
phase for this model 0, [|[ , it has been shown that the 
liquid pocket in the Cgo phase diagram is confined to a 
tiny temperature interval (see e.g. @, 0] and references 
therein). In this sense, the system displays a characteris- 
tic "borderline" behavior, intermediate between what one 
expects for the phase equilibria of a simple fluid (with a 
fully developed liquid phase), and a condition where the 
liquid-vapor equilibrium is only metastable with respect 
to the vapor-solid phase separation, the binodal curve 
falling below the sublimation line. The latter behavior 
is usually observed when the range of attractive forces is 
short enough compared with the size of the repulsive core, 
a condition typically faced when one considers effective 
models for macrosized molecular systems, like protein so- 
lutions d, HI, |T(| or colloidal suspensions [ill • 

On the other hand, extensive studies of the glass tran- 
sition in simple systems like square wells, adhesive hard 
spheres and hard-core Yukawa fluids [lj, 03 El ( see 
ref [l5[ for a recent review), based on the application 
of the Mode-Coupling Theory (MCT, 0) have allowed 



to identify distinct mechanism of formation of glasses 
according to the balance of the repulsive and attrac- 
tive interactions. In particular, systems with sufficiently 
short-range attraction exhibit, together with a normal 
repulsion-driven glass which behaves qualitatively like 
a hard-sphere glass, an "attractive" glass of different 
nature, favored both by the energy and the local en- 
tropy [13, EH- This circumstance naturally candidates 
the C 6 o model — wherein, as observed, the subtleties re- 
lated to the shape of the interaction potential play a cru- 
cial role — for a study of the glass transition and of its 
typicality. 

Besides our simulations investigations @, [|[ , the study 
of the glass transition line in the Girifalco model has been 
recently addressed by Greenall and Voigtmann In 
ref [l!| these authors carry out ideal MCT calculations, 
using as input data for the theory the static structure fac- 
tors S(k) obtained from the Hypernetted Chain (HNC) 
and Percus-Yevick (PY) liquid state theories [2(j. They 
show that vitrification in the Cgo model occurs, in agree- 
ment with MD results, although at densities lower than 
those predicted by the computer simulations; such an 
underestimate is not unexpected given the inherent in- 
accuracies of MCT [14]. Moreover, the features of the 
MCT non-ergodicity parameter and the overall behavior 
of the glass transition line indicate the crossover to an 
attractive-glass behavior at relatively low temperatures, 
thereby expanding the scenario emerging from previous 
simulations. MD simulations have shown in fact evi- 
dences of a repulsive glass only, over the whole temper- 
ature range investigated 0, Q ■ The possibility that an 
attractive glass can exist for the Girifalco model appears 
somewhat unexpected on the basis of the broad analysis 
carried out in ref I14j ]. in which a Yukawa model with 
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parameters compatible with the decay rate of the Giri- 
falco potential displays a repulsive glass only. A similar 
conclusion is drawn in ref fl9 ] itself: if one uses a square- 
well potential mimicking the attraction range of the Gir- 
ifalco model, the resulting interaction is not short-range 
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enough to determine the appearance of an attraction- 
driven glass. 

The possibility that the attractive interaction in the 
Girifalco model is sufficiently short-range to display sev- 
eral peculiar features of an attractive glass poses intrigu- 
ing questions about the whole mechanism which under- 
lies the existence of such an arrested state of matter, 
with implications for a variety of similar potentials cur- 
rently used in colloid and protein studies. On the other 
hand, the arguments put forward in ref [l9| hinge on two 
approximate liquid state approaches, like the HNC and 
the PY theories, and on a previous analysis carried out 
by the same authors [2l|, about the (weak) sensitivity 
of MCT predictions to the S(k) behaviour in the low- 
k region, below the first diffraction peak. This scenario 
motivates in our opinion an investigation of MCT pre- 
dictions implemented through the use of accurate struc- 
ture factors. We have performed to this aim extensive 
MD calculations in the temperature and density regimes 
inside the vapor-solid and liquid-solid regions, and ex- 
tracted the structure factor directly from the simulation 
data, down to the lowest k vector compatible with the 
simulation box size. In order to assess the theoretical 
predictions, we have also calculated the static structure 
factor in some selected thermodynamic state points, in 
the framework of the the Modified HNC (MHNC, J^) 
approach, a theory which yields accurate results for the 
C 60 model (H. 

The paper is organized as follows: sect. II is devoted 
to an introduction of the model, of the simulations strat- 
egy and of the basic equations employed in the MCT. 
In sec. Ill results are reported and discussed, whereas 
sect. IV contains a few concluding remarks. 



II. MODEL, SIMULATION PROCEDURE AND 
MCT APPROACH 



The Girifalco interaction potential v(r) among two Cgo 
molecules reads nf: 
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where s = r/d, on = N 2 A/I2d 6 , and a 2 = N 2 B/90d 12 ; 
N = 60 and d — 0.71 nm are the number of carbon 
atoms and the diameter, respectively, of the fullerene 
particles; A = 32 x 10~ 60 ergcm 6 and B = 55.77 x 
10~ 105 erg cm 12 are constants entering the 12-6 potential 
</>(r)= — A/r 6 + B/r 12 through which two carbon sites 
on different spherical molecules are assumed to inter- 
act. The distance where the potential in eq [1] crosses 
zero, the position of the potential well minimum and 
its depth, are a = 0.959 nm, r m - ln = 1.005 nm, and 
e = 0.444 x 10~ 12 erg, respectively. 



MD simulations are carried out on a system composed 
of 1000 particles enclosed in a cubic box with periodic 
boundary conditions. The Andersen algorithm 24] is 
used to enforce constant-pressure P, constant-enthalpy 
H conditions to the sample. We have analyzed the sys- 
tem along the cooling paths characterized by the pressure 
P = 3,5, 40, 150 and 250 MPA^iccording to the strategy 
already documented in refs 2, 3] . As visible in Figure [TJ 
five to seven thermodynamic states for each pressure are 
selected around the glass transition points (which have 
been estimated through several structural and dynamic 
indicators in ref [3J]), and long trajectories are generated 
to estimate the equilibrium static structure factor S(k) to 
be fed into the mode-coupling theory calculations. Sev- 
eral state points along the isotherms T — 1200 K and 
T = 3500 K are also analyzed, in order to characterize 
the behavior of the vitrification curve in the limits of 
relatively low and high temperatures (see Figure [T]) . 

Runs of 12 000 time steps (with At = 5fs) are gener- 
ally long enough to ensure a stable estimate of structure 
factors. Two method for the calculation of S(k) (namely 
through a direct estimate of fluctuations of the density p 
and by fourier inversion of the radial distribution function 
g(r)) are used and compared in this study. Results are 
generally equivalent: the calculations coming from g(r) 
are on the whole less noisy than those obtained through 
the direct method whereas the latter are more accurate in 
the small-A; region. If necessary, a smoothing procedure 
has been applied to the S(k) data prior to mode-coupling 
calculations. 

MCT derives equations for the normalized time- 
dependent density correlators of the Fourier components 
of the particle density fluctuations Spk(t) 16 : 



(6p k ®5p- k (0)) 
\5p k 
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starting only from the number density p and the struc- 
ture factor Sk = (\Spk\ 2 ) /N (or equivalently the direct 
correlation function pek = 1 — S^ 1 . The glass transition 
predicted by MCT is obtained solving the t — > oo limit 
of the equations for the normalized correlators $k (t) , the 
so-called non-ergodicity factor 
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The equations have the form: 
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FIG. 1: Cooling paths at P = 3.5, 40, 150, and 250 MPa 
and constant temperature runs at T = 1200 and 3500 K. Full 
lines correspond to the liquid branch of the binodal and the 
liquid-solid coexistence lines [6J. The triple point (diamond) 
and the metastable portion of the binodal line (dashed curve) 
are also shown. Dots collectively show all past [(| and newly 
added (see text) estimates of coexistence points. 



For specific values of the input parameters, the solution 
to these equations admits not only the usual trivial so- 
lution fk = 0, but also solutions with fk ^ 0. The value 
of fk at the transition point is denoted /£. We have 
solved equations [4][5] on a grid of 300 wavevectors up to 
k = 32 nm -1 using a standard iterative procedure. 



III. RESULTS 

An overview of all state points encompassed in this 
work is reported in Figure [TJ in the context of the phase 
diagram of the Cgo model calculated in ref Newly 
generated fluid-solid coexistence points in the high tem- 
perature regime (T — 3500) and in the low temperature 
solid phase (T < 1200) are displayed as well, in order to 
elucidate the whole appearance of the coexistence region 
where all calculations have been done. The coexistence 
points at high temperature are in particular calculated 
according to the procedure employed in ref (2f| : the free 
energy of the fluid phase is calculated through thermo- 
dynamic integration of the MD pressure at several state 
points along the isotherm T = 3500 K, whereas for the 
solid phase we have used a first-order perturbation the- 
ory starting from a crystal of hard spheres, whose diame- 
ter is chosen accor ding to the Weeks- Chandler- Andersen 
(WC A) procedure [26[ . We have obtained in this way for 
the coexisting densities: pfl U id(T = 3500) = 1.02 nm -3 
and p so \id(T — 3500) = 1.24nm -3 . As for the isotherms 
T = 1200 K and T = 800 K, wc have assumed that the co- 
existing solid density must be that where the (perturba- 
tion theory) pressure reduces to zero — a simple approxi- 
mation justified by the fact that in this regime the coex- 
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FIG. 2: MD static structure factors along the T = 1200 K 
(left) and P = 3.5 MPa (right) paths. In the top panel, the 
corresponding non-ergodicity factors are also displayed. 



isting vapor density is extremely low and thus it plausibly 
corresponds to almost zero pressure. We have obtained: 
p so m{T = 1200) = 1.34 nm- 3 and p so iid(T = 800) = 
1.36 nm -3 . 

As visible from Figure [TJ thermodynamic states at 
pressures increasingly higher than P = 3.5 MPa are allo- 
cated between the freezing and the melting lines of the 
model. State points along the isobaric path P = 3.5 MPa 
are almost superimposed on the liquid branch of the 
metastable portion of the liquid-vapor phase separation, 
whereas all thermodynamic states along the isotherm 
T = 1200 K are definitely inside the vapor-solid coex- 
istence region. This point can be further illustrated by 
the evaluation of the structure factors at T = 1200 K 
and P = 3.5 MPa, reported in Figured It appears that 
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FIG. 3: Snapshot of the final configuration of the system at 
T = 1200 K and p = 0.95 nm -3 ; a void region clearly displays 
in the sample. 




time [ns] 

FIG. 4: Mean square displacements across the MCT vitri- 
fication thresholds at P = 3.5 and 250 MPa. Longer runs 
correspond to the states immediately before the transition 
point. 



the structure factor at T = 1200 K is barely sensitive to 
density variations in the range 0.95 — 1.05 nm -3 , but for 
the steep rise of the k — > limit when the density is de- 
creased, clearly indicating a more and more pronounced 
tendency of the sample to phase separate as it goes deeper 
and deeper inside the (metastable) binodal. On the other 
hand, the structure factor at P = 3.5 MPa keeps around 
zero in the k — > limit and becomes more structured 
during the cooling sequence documenting in this way the 
approach to the solid configuration. 

In parallel, the spatial configurations obtained at T = 
1200 K, analyzed in terms of nearest-neighbors distribu- 



tions of particles and by means of a direct visual inspec- 
tion, reveals at p = 0.95 nm -3 the existence of a rar- 
efact region in the sample surrounded by a uniformly 
denser environment (see the snapshot in Figure [3]) . Such 
strong dishomogeneities tend to persist almost up to 
p = 1.025 nm -3 and disappear at higher densities: at 
p = 1.05 nm -3 we have detected no evidence of (incipi- 
ent) phase separation, and the system displays a uniform 
distribution of particles inside the simulation box. 

At variance with the latter regime, all other points in- 
vestigated in this work represent truly equilibrated ther- 
modynamic states. In order to be sure that all sam- 
ples are properly equilibrated, and that each particle dif- 
fuses on average more than its diameter over the whole 
simulation time, we have recorded the mean square dis- 
placement (MSD) over particularly long simulation runs 
for some selected state points where MCT predicts the 
glass transition. Results are reported in Figure 2) as vis- 
ible, observation windows of the order of ~ 0.1ns allow 
the particles to diffuse for a distance varying between 
~ 1.4 and ~ 1.7 molecular diameters, the latter being 
roughly identified with the "collisional" distance a of the 
Ceo potential (see section II). The diffusion coefficient 
varies between ~ 2.6 x 10~ 5 and 3.8 x 10~ 5 cm 2 /s be- 
fore the MCT transition and drops to ~ 1.6 x 10~ 5 - 
6.6 x 10 _6 cm 2 /s thereafter. Data in Figure S] confirm 
the known MCT tendency to predict a glass transition 
at state points where the dynamics of the real system is 
still ergodic. 

Main results of this work are presented in Figure [5l 
where the glass transition line identified through MCT 
calculations is displayed. We also report in the same fig- 
ure the theoretical predictions from MHNC calculations 
at T = 1200 K and 3500 K, our previous MD data for vit- 
rification @,[3j|, and the estimate of ref [H|. MHNC cal- 
culations follow the general approach already employed 
successfully to predict the liquid-vapor coexistence of the 
Ceo model at issue 23]. The numerical solutions of the 
MCT long-time limit equations (the non-ergodicity fac- 
tors) across the vitrification thresholds are shown in Fig- 
ure [6l for several isobaric paths investigated in this work. 
The last ergodic state points (//. = 0) and the first non- 
ergodic ones (/& ^ 0) are assumed to bracket the ideal 
MCT glass transition line reported in Figure The dis- 
tances in the p-T plane among these couples of points 
hence constitute the error bars of our predictions. 

It appears from FigureOthat the MCT glass transition 
line does not show re-entrant behavior and is reasonably 
parallel to the vitrification locus obtained in our previ- 
ous MD simulation study 0, [1[ (we shall further com- 
ment about these specific features of our results), and 
runs fairly close to the one determined through MCT and 
HNC or PY structure factors as input data in ref [l!| . In 
particular, PY predictions tend to slightly underestimate 
the vitrification density and the HNC results are gener- 
ally more accurate over the temperature regime investi- 
gated. 

Our test MHNC and HNC calculations, limited to 
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FIG. 5: Glass transition points as obtained in this work 
through MCT calculations based on MD (full squares), and 
MHNC (triangles) structure factors. Crosses are the glass 
transition estimates through direct MD calculations [3, 01 • 
MCT predictions with HNC (circles) and PY (open squares) 
structure factors [19| are also shown. 
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FIG. 6: Non-zero ergodicity factors at the MCT transition 
thresholds (/£) for pressures P = 250, 40 and 3.5 MPa (from 
top to bottom). State points along such isobaric paths with 
temperatures immediately higher than those reported in the 
figure are still ergodic, i.e. fk = 0. Densities and tempera- 
tures in the legends are expressed in nm -3 and K, respectively. 



some selected densities along the isotherms T = 3500 K 
and T — 1200 K are also reported in Figure [5] and con- 
firm the trend illustrated above. A comparison between 
theoretical and simulation structure factors is reported 
in Figure El We display in particular data at the high- 
est temperature and density investigated in this work, 
namely T = 3500K and p = 1.15nm -3 , and data at the 
opposite extremum, i.e. T = 1200 K and p = 1.05 nm -3 
(the lowest state point which we have been able to inves- 



tigate before the system displays a pronounced tendency 
to phase separate). Apart from the well-known difficul- 
ties of the theoretical tools employed to closely follow the 
k — > behavior of the structure factor, especially in the 
proximity of the phase separation, it appears that both 
theories slightly underestimate the height of the main 
peak (with MHNC theory performing better than HNC) , 
and a small dephasing emerges in the position of the sec- 
ondary peaks at higher k vectors. On the whole, we judge 
the agreement between theory and simulations quite sat- 
isfactory, especially if we take into account that all data 
recorded in Figure [7] have been obtained by pushing the 
integral equation scheme toward "extreme conditions" , 
as far as the regime in which such tools are known to 
give the best performances is concerned. 

Figures[5]and[7]document the good agreement between 
present simulation results, our thoretical predictions, and 
those obtained by Greenall and Voigtmann equally 
from the point of view of MCT transition threshold, and 
at the level of input structure factor calculations. Our 
evidences indirectly support the conclusions drawn in 
ref [l9| that, in the framework provided by the mode- 
coupling theory, and in the thermodynamic region T < 
3000 K, the attractive part of the Girifalco potential is 
short-range enough to begin to influence the approach 
to the structural arrest of the system. The arguments 
presented in ref [l|| are based in particular to an ex- 
tended analysis of the properties of the non-ergodicity 
factor, together with a comparison between the MCT 
lines obtained with either the full Cgo potential, or with 
a truncated version, where the attractive part is cut off 
according to the WCA procedure. 

We now turn to the examination of the low- 
temperature regime MCT results. As visible in the top 
panel of Figure at T = 1200 K MCT predicts the ex- 
istence of a glassy phase for p — 1.05 nm~ 3 , which cor- 
responds, as discussed above, to the lowest density for 
which a genuine homogeneous sample is observed during 
the time of our MD simulations, whereas the tendency 
of the system to phase separate is already well under 
way at p = 1.025 nm -3 . For this reason, and in view 
of the fact that MCT itself strictly deals with homoge- 
neous systems, we estimate the vitrification density at 
T = 1200 K as falling halfway between p — 1.025 nm -3 
and p = 1.05 nm~ 3 ; this value, reported in Figure[5]nicely 
coincides with our MHNC estimates. If lower density — 
but phase-separated — MD structure factors are fed into 
the MCT equations, glass states are equally predicted. 
The reason of such an outcome can be appreciated if we 
observe, again in Figure^ the way the differences in the 
shape of the various structure factors reflect in the corre- 
sponding non-ergodicity factors: it appears that, in pass- 
ing from p — 1.050 nm -3 to p = 0.950nm~ 3 , the struc- 
ture factor displays minor differences, (especially in the 
sensitive region of the nearest-neighbor peak) , except for 
the already commented steep rise at k = 0. By contrast, 
the non-ergodicity factors, also shown in the same figure, 
are hardly affected by the density variation, and infact 
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FIG. 7: Comparison between theoretical and simulation re- 
sults for the static structure factors at T — 3500 K and 
p = 1.15 nm -3 (top) and T = 1200 K and p = 1.050 nm" 3 
(bottom). 



run almost superimposed on top of each other, thus pre- 
dicting a glassy phase all along the isotherm T — 1200 K. 
These evidences support the theoretical observation that 
MCT predictions for glasses characterized by short-range 
attractions, are barely sensitive to the low-fc behavior of 
the structure factor [2X1 ] . 

The non-zero ergodicity factors displayed in Figure [2] 
would allow us to shift to lower densities the T = 1200 K 
MCT vitrification threshold, causing a net bending of 
the glass transition line, even more pronounced than that 
predicted in ref (l9( . Although this feature would indicate 
an enhancement of the attraction-driven properties of the 
glass phase, the tendency to phase separate of our sam- 
ple detected during MD simulations prompts us to some 
caution as far as any definite conclusion about this point 
is concerned. Similar problems with the phase separation 
of the sample prevented us from investigating the vitri- 
fication threshold for temperatures T < 1200 K. More 



TABLE I: Comparison, for several temperatures, among the 
MCT transition densities calculated in this work (column 2), 
the MD results of refs |2j,[3J (column 3) and the estimate which 
are obtained if we apply a 12% shift to MCT data (colunm 
4). Temperatures are given in K, densities in nm -3 . 



T 


MCT 


MD 


MCT (shifted) 


1635 


1.073 


1.190 


1.206 


1478 


1.050 


1.187 


1.180 


1200 


1.037 


1.177 


1.166 



extensive simulations, accompanied by an analysis of the 
size dependence of the results, might be worth in this 
respect. Such a program however is unlikely to be imple- 
mented at ease, given the fact that the thermodynamic 
region to be explored falls well beneath the metastable 
liquid- vapor binodal line of Cgo, and definitely inside the 
solid-vapor metastable equilibrium region, both circum- 
stances implying a strong trend to phase separation. 

In comparison with our previous MD estimates of the 
glass transition of the Girifalco Ceo model based on 

several structural and dynamic indicators, the MCT pre- 
dictions moderately underestimate the transition densi- 
ties. It is well known that a corresponding underestimate 
occurs when reference MCT predictions for hard spheres 
(%iass = 7r/6/5fj 3 = 0.516, [271 ]) are compared with the 
rigorous computer simulation result (?7giass = 0.58, [281] ) . 
In this case, a common procedure has been to exploit 
such a ~ 12% discrepancy, in order to sistematically 
shift to higher densities the MCT vitrification curves, 
also for other models, when a comparison with simula- 
tion data is carried out (see e.g. ref [HI). If we ap- 
ply the same correction to our MCT points, we obtain 
the remarkable agreement with MD estimates reported 
in Table Q] In refs 0, 0| we have shown that all vit- 
rification densities reported in Figure [5] nicely fall on a 
single value ((Tr/6pa^ s ~ 0.58, almost coincident with 
the hard-sphere glass transition density), on condition 
that all data are rescaled through the effective (tempera- 
ture dependent) diameter a e s which is obtained by sub- 
stituting the soft repulsive part of the Girifalco poten- 
tial with a pure hard-core exclusion, according to the 
WCA prescription. As a consequence, we have deduced 
in refs [3, Q that the hard-sphere behavior plays the main 
role in determining the structure of the Cgo glass, which 
hence results essentially repulsion-driven in nature. In 
order to reconcile the attractive nature of the MCT glass 
discussed in ref with the results of refs 0, Q it is 
worth observing in Figure [5] that in the low tempera- 
tures regime the bending of the MD curve is milder, with 
respect to the MCT predictions. In this case, it could 
be possible that, on lowering the temperature further, 
also the glass observed in our previous studies would dis- 
play a more pronounced attraction-driven effects. On the 
other hand, the analysis carried out in this work indicates 
that the strong tendency of the system to phase separate 
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would likely preclude the possibility of extracting mean- 
ingful and robust conclusion from an investigation of the 
behaviour of the Ceo model in this (highly metastable) 
low-temperature, high-density regime. 

IV. CONCLUSIONS 

We have reported MCT determination of the glass 
transition line in a model Cgo- The theoretical cal- 
culations are based on the use, as input data for the 
MCT equations, of structure factors obtained via exten- 
sive molecular dynamics simulation investigation of the 
fullerene systems at various pressures. The shape of the 
vitrification locus is in fairly good agreement with our 
previous determinations of the glass transition entirely 
based on MD simulations @, In fact, the MCT pre- 
dictions appear to overall underestimate the glass den- 
sity by ~ 10% with respect to refs |2j, |3|, manifesting an 
otherwise well-known inaccuracy of the theory. Present 
results are also in fairly good agreement with the glass 
transition line obtained by other authors through MCT 
calculations based on HNC structure factors input [l9j |. 
MCT calculations based on the refined modified HNC 
theory S(k), also presented in this work, also reproduce 
faithfully the glass transition line. 

As far as the interplay between the non-ergodicity fac- 
tors and the structure factors S(k) is concerned, our evi- 
dence lends support to the theoretical analysis of ref , 
where the influence of the attractive part of the Ceo po- 
tential on the structural arrest properties of the model is 
argued. Our simulations document a strong tendency of 



the system to phase separate in the low-temperature re- 
gion where the MCT glass transition line should exhibit 
a more pronounced attraction-driven character. In this 
case, as required by MCT, we have restricted our analy- 
sis only to those temperatures and densities for which a 
fully homogeneous sample is obtained, this choice result- 
ing in a glass line with rather a moderate bending. More 
extensive computer simulation investigation, which take 
into account for instance size effects, might be worth in 
order to further enlight this point. 

The picture emerging from the overall shape and 
characters of the glass transition line vs the repulsive- 
attractive potential features turn out to be consistent 
with the more general proposition that the glass line 
terminates on the liquid side of the liquid-gas coexis- 
tence, when particles interact with spherical potentials 
in which the excluded volume repulsion is complemented 
by attraction [29|, [3(|. This scenario actually appears 
to be valid independently from the range of the attrac- 
tive potential [291 ]. Non-spherical patchy potentials, in 
which the number of interacting particles is significantly 
reduced as compared to the spherical case, are neces- 
sary in order to suppress and shift to small densities the 
liquid-gas phase separation curve [3ll . I32I ] and hence to 
extend the dynamic arrest line to lower temperatures and 
smaller densities. 
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